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Abstract 

Inelastic lifetime of an electron quasiparticle in an electron liquid due to electron-electron in- 
teraction evaluated in previous work is calculated in an alternative way. Both the contributions 
of the "direct" and "exchange" processes are included. The results turn out to be exactly the 
same as those obtained previously, and hence confirm the latter and consequently fully resolve the 
theoretical discrepancies existing in the literature. Derivation in the two-dimensional case is pre- 
sented in detail due to its intricacies. The effects of local field and finite well width on the effective 
electron interaction in the two-dimensional case are also investigated in a quantitive comparison 
of the electron relaxation rate between theory and experiment. These effects are shown to make 
rather small contribution to the quasiparticle lifetime. 

PACS numbers: 71.10.Ay, 71.10.-w, 72.10.-d 
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I. INTRODUCTION 



Low-energy electron excitations in solids can be successfully described in terms of quasi- 
particles in the Landau theory of Fermi liquids. 

ELQ 

An excited quasiparticle with definite 
momentum is not stable due to scattering by phonons, disorders, and other electrons. Hence 
a quasiparticle has finite lifetime. Among this, the intrinsic inelastic scattering lifetime r e , 
i.e., the lifetime that arises purely from the electron-electron scattering processes, is a central 
quantity in the Laudau theory of the electron liquid. It plays a key role in our understanding 
of a broad variety of phenomena in solids such as electron dephasing, (s| tunneling, jj, 0] and 
localization,!™ etc. It might also have effect on electron transport. 17] 

nn 

In fact, the electron tunneling techniques in semiconductor quantum wells |5( have 
enabled experimentalists to directly determine r e in two-dimensional (2D) electron liquids. 
For weakly coupled wells, the lifetime principally arises from electron-electron scattering 
processes. On the other side, huge progress has also been made, by the use of the techniques 
of ultrafast laser, in measuring the lifetime of photonexcited electrons in metals such as 
copper. P These advances have made it possible to carry out quantitive comparisons be- 
tween theories and experiments. The theory of the inelastic lifetime in three dimensions 
(3D) is rather well established within random phase approximation. P. llo| It was later ex 
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tended to include the exchange contribution (see Ref. ll| for a detailed discussion). 
Several theoretical investigations had also been carried out in 2D, but with quantitive 
disagreement. In an earlier paper, [111 hereafter referred as I, 

we have managed to clarify the origin of the disagreement that exists among these previous 
investigations. The results in I are summarized as follows. The inverse lifetime of a quasi- 
particle with low energy £ p (relative to the chemical potential /z) at temperature T in a 3D 
electron liquid is 



1 m 3 e 4 v 2 k 2 B T 2 + e p 



r e irpkf 1 + e ftp 



A 



+ tan 1 A 



A 2 + 1 

(1) 



7T _J\ I 1 

— tan 



y/J? + 2{2 VAVA 2 + 2 

where A = 2kp/k s , and k s = yj is the 3D Thomas-Fermi screening wavevector. ks, kp, 
and ao are the Boltzmann constant, the Fermi wavevector, and the Bohr radius, respectively. 
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For a 2D electron liquid, we found 



1 m 2 £ 2 

^ p r 3W 2 {0) + 2W 2 {2k F ) 



r e 167r 3 £ F ' 



for fcgT <C £ p , and 



j_ (mk B T) 2 
T e ~ 32nE F 



-2W(0)W(2k F )]ln^r (2) 



3H/ 2 (0) + 2H/ 2 (2A; F ) 

-2M/(0)H/(2A; F )]ln^ (3) 



for £ p <C A;_bT. Here Ep = h 2 k F /2m, and W(q) is the effective interaction between quasi- 
particles. 

The calculation of r e is a quite nontrivial task in many-body theory, which helps explain 
the disagreement among various previous theoretical results in 2D. Evidently, correctness of 
the results in Eq. (P), and Eqs. (J2J) and (J3j) is crucial in any meaningful comparisons with 
experiments. In fact, big dis crep ancies remain between experiments and theories, and call 



tact, big discrepancies remain between 



for explanations.^, y, la, lUl, [lj, [la, [ly, [TO, llSl llU l2l| In this paper, we calculate l/r e in 
an alternative way, and confirm the above results. Various theoretical discrepancies in the 
literature, as mentioned above, are fully resolved. 

The present calculation appears very different to the previous calculation in I. In I, we 
calculated l/r e by expressing it as the frequency convolution of the imaginary part of the 
density-density response function. The present calculation is technically straightforward, 
and somehow in a textbook fashion. But it is by no means much simpler than that in I. In 
fact, the present approach seems rather clumsy for the case of £ p k B T ' . Hence we shall 
restrict ourselves to the case of zero temperature. 

After giving the general formulas for l/r e in the next section, we present our calculation 
for the 3D and 2D cases separately in Sec. Ill and IV. In Sec. V, we shall discuss the 
contribution to the inverse quasiparticle lifetime in 2D arising from the effects of local field 
and finite well width on the effective electron interaction, and then briefly summarize the 
paper. 
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II. GENERAL FORMULAS 



We start by rewriting Eqs. (4) and (5) in I as follows: 

-J5j = 271 W ™' ( k - P)™ka™k+q-pa'™qa' 

T °" k,q a' 

^ (f p + f k+q-pa' - £ka ~ £q CT ' ) > ( 4 ) 

and 

= -2vr ^Up " q)W w (k - p) 

k,q 

^k<7^k+q-pcr''^qcr'<K£p + £k+q-p<7 — £k(j ~ £q<r) i (5) 



where W aa r(q) is the effective interaction between quasiparticles of spin a and spin a'. We 
have set H — 1. We consider only the paramagnetic electron liquid, and hence the index a of 
the l/Ta^'^ is unnecessary and will be dropped hereafter. The sum of and I/t^ 

yields the inverse inelastic lifetime l/r e , 

- = 27r EE (l-^)[W^(k-p) 



k,q a' 



Saa'Ww'ip - q)] 2 n k(7 n k+q _ p(7 



& (f p + 6c+q-pa' - f ~ W' ) ■ ( 6 ) 

We are only interested in the case of low excited energy (£ p <C Ep). In this case the 
contribution to the summations over momenta in the above expression arises only from the 
region of £ k +q- P <7, £k CT , £qa < Therefore 

— = 2tt ^2^1 r<J ' / i q) r2 ka^k+q- PCT 'n q(7 ' 

^ k,q a' 

£(£p<7 + £k+q-p CT ' - £k CT - £qa')> (?) 

where /x k — p-k and /i q = p-q respectively, and the hats mean unit vectors. We have defined 
formally 

-S^W^iVlkF^/T^y)] 2 . (8) 

Below we present our calculation for the 3D and 2D cases separately in Sec. Ill and Sec. 
IV. We set the volume of the 3D system and the area of the 2D system, respectively, to be 
unit in this paper. 



III. THE INVERSE LIFETIME IN 3D 



The integrations over the azimuthal angles of k and q with respect to p in Eq. (J7j) can 
be straightforwardly carried out. After that, it becomes 



9(p 2 + kp - k 2 - q 2 ) J dfik J d/i q 



where £ p = p 2 /2m. We have defined B(x) = 1 far * > 0, 0( X ) = for * < 0. Notice that 
Eq. Q is essentially the same as Eq. (4) in Ref. |l3j] except typos of (^sO 3 and a missing 
factor of 1 — \5 aa i in Eq. (3b) in Ref. [13]. The integrations over k and q in Eq. (j^J) can be 
carried out analytically to the accuracy of the leading order of 0(^ 2 ), and yield 

i . £ 



where 



^(,- L ,,-J , 1 ^ + ^) (io) 

q V(i-^)(i-^ q ) 



fl 



Notice that there seems an error of a factor 1/2 in Eq. (6b) in Ref. However this 

does not effect the results in Table II in Ref. |13J] and the subsequent conculsions, since the 
coefficient A cancels in the ratio p aa /p cra in the table. 

The contributions from the "direct" and "exchange" processes are, separately, given as, 



1 £2 rl rl 

— — - A ' 

r( fl ) E F E 



[H/K^v^T^) + W*,(y/2k P y/T=]£)] 
1 1 



V"(l - /ik)(l - /iq) V^k + /V 



(12) 



where cr = — a, and 



1 £ 2 f 1 f 1 

1 1 
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(13) 



The spin dependence of the effective interaction is not crucial in determining the to- 
tal inelastic lifetime. |3] We shall ignore this dependence and follow the usual practice of 
characterizing the screening effects by the screening wavevector k s : 

4ire 2 



W(q) 



(14) 



q 2 + k 2 

The integrations over /x q and /Xk can be analytically carried out. After that, one finally has 



m 3 e A k 



F t 2 



T 



(D) 



A 



A 2 + l 



+ tan 1 A 



(15) 



and 



m 3 e A k 



T 



(ex) 



np 2 k 3 s " p ^W+2 



- - tan" 1 ( - 



1 



A V A 2 + 2 



(16) 



Equations (|T3|l and (fTB^l are, to the leading order of 0(£ 2 ), exactly the same as the results 
for 1 / t( d ) and 1 / r^ ex ' obtained in I for the 3D case. 



IV. THE INVERSE LIFETIME IN 2D 

The derivation in the 2D case is relatively intricate. We shall present it in details. We 
first rewrite Eq. (J7J) as 

1 / 1 \ 2 f°° f°° 

— = 2irm ( - — — - I > / dqq / dkk 
\(2n) 2 J ^J kF qq J kF 

/TV r>TT 
dfa / # q I™'Gu k ,// q ) 
-7T J —7T 

5(p 2 + kq cos(0k — q ) — pk cos 0k — pq cos q ), (17) 

where /ik = cos0k, /^q — cos0 q , and 0k and q are the angles of k and q relative to p, 
respectively. After some algebraic manipulations, the preceding expression can be written 

as 



1 / 1 \ 2 TP ry/p 2 +k 2 F -q 2 

J dfi k J rf/i q / CTCT '(/i k ,Ai q ) 

6[(p ~ k^) 2 (p - q^f - k 2 q 2 (l - /4)(1 - £)]. 



T, 



(18) 



y+(x) A 2 

A 



D 



FIG. 1: The contour of f(x,y) = 0, denned in Eq. 1)20(1 . and its two parts: y+{x) (solid line) and 
y_(x) (dotted line), defined in Eq. (|28|). 



We now define the following variables: 



and function: 



(19) 



f{x, y) = (A - x) 2 (A' - j/) 2 - 1 + x 2 + y 2 - xV 



(20) 



and write Eq. (|18|) in the following form: 



1 2-/// / 1 

r P . 



/t 2 Ik 2 

cr' K F ° K F 



dz' 



dx [ dy6[f(x,y)]r^(x,y). 
-l J-i 



The contour of f(x,y) = is illustrated in Fig. 1. 
We further define 

CP 



siA,B,D 



-p 2 +k 2 F -z 

dz 



k'p J k^ 



7 / jA,B,D 



where 



J ao' 



dx [ dy6if(x,y)]r°\x,y), 
o Jo 



(21) 



(22) 



(23) 



J*,= / dx I dy6[f(x,y)]I™'(x,y), 
1 Jo 



/•! 



(24) 
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and 




dx / dy5[f(x,y)]I^(x,y) 



-i 



(25) 



and express l/r e as follows: 



1 



"^((2^) E^-' + G -' + G -'l 



(26) 



Evidently, A, B and D denote the contributions arising from the first, the second and the 
fourth quadrants respectively, as shown in Fig. 1. Notice that f(x,y) ^ in the third 
quadrant. In fact, the leading order contributions only arise from the circled regions in Fig. 
1. By interchanging the integral variables x and y in Eq. (j25|) . it is easy to see that the 
evaluation of G® a , is totally analogous to that of G^ a ,. Hence we only present the latter. 
The evaluation of G^ a , turns out to be relatively simpler than that of G^ a , , and it is in the 
meanwhile instructive for the latter. We hence start with the former below. 

We first evaluate the integrations over the variables x and y in Eq. (|24j). To this end, we 
rewrite Eq. (J2"lj) as 



The functions y+{x) and y~(x) are two components of the contour of f(x,y) = 0, and 
they are illustrated in Fig. 1. The integration over the variable y in Eq. ()27|) is now 
straightforward, which yields 




6([y-y + {x))[y-y4x)})I^{x,y) 



(27) 



where 




(28) 




2 7-1 v / T^^v / ( 1 - A ' 2 )( A -x) 2 + l-x 2 
[I^(x iy+ (x)) + I^(x iy 4x))}, 



(29) 



or, more explicitly, 




[I^(x,y + (x)) + I^(x,y4x))}. 



(30) 
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Here we have defined 



A(l - A /2 ) ± VA 2 + A' 2 - A 2 A' 2 
^1,2 = y 2 , (31) 



which are also shown in Fig. 1. By using the fact that xi — 4(A' — 1) — 1 and y + (xi) = 
U-(xi) = 1 for A' — > 1, we have, to the leading order, 

J^, = -V'(-l,l)ln[A'-l]. (32) 

Substituting the preceding result into Eq. (|22j) and performing the integrations over z and 
z', one obtains 



Gt> = -\ iaa '(-^ !)(P 2 " k 2 F ) 2 H(P 2 ~ 4)/24l- (33) 

As pointed out previously, the evaluation of G® a , is similar to that of G^ a ,. Here we only 
quote the final result, 

= -1)(P 2 ~ k%f ln[( P 2 - 4)/2fc 2 ]. (34) 

Next we calculate G^ a ,. First of all, from the experience in deriving J^ a , in Eq. ()32|). it 
is not difficult to see that the leading order contribution to J^ a , arises from the region of 
x — > 1, y — > 1. Therefore, we may directly rewrite Eq. (|2"3*|) as 

J CT A CT , = I™' (1,1) [ dx [ dy5[f(x,y)}. (35) 
Jo Jo 

However, the following calculation is a little more delicate. Due to the fact that y~(x) be- 
comes ill-defined (actually becomes x — 1) in the first quadrant as A, A' — > 1, a straightfor- 
ward calculation like the preceding one for J^ a , does not work. To circumvent this difficulty, 
we make the following variable transform, 



X = l=( X '-yi), y =lA.( X < + y>). (36) 



and rewrite Eq. (J3l>j) as 



j£,= 7^(1,1) y dy' 

8{a{x')[y'- yi {x')][y'-y2{x')]). (37) 

The Jacobian A'/A for the above integration variable transform can be set to be one in the 
limit of A, A' — > 1. The limits of the integration over y' are left unspecified because they 



are not really relevant simply due to the ^-function in the integrand, while the integration 
region of x' and y' corresponds to the square of < x', y' < 1. In Eq. (j3*Tj) . we have defined 



where 



-b(x) ± ygg) - 4a(x)c(x) 

^ 2(x) = Mx) ' (38) 



a{x) = ^[2V2x\ a /\ - 2A' 2 + A /2 /A 2 + 1], (39) 



b( x ) = -(1 - A ,2 /A 2 )x, (40) 



and 



c(x) = A 2 A' 2 - 1 - 2v / 2AA ,2 x 
1 



+-(A ,2 /A 2 + 6A' 2 + l)x 2 - V2(A /2 /A)x 3 . (41) 
2 



The integration over y' yields 



Jo v 6 2 (x) — 4a(x)c(a;) 



which can be rewritten as 



(42) 



^=/-(M)/° *y, + * +7) . (43) 

7--J=(i+a/a') Var + + 7 
where a, /3, and 7, in the limit of A — > 1, A' — > 1, can be shown as 

a = 16, (44) 

(3 = -4y/2[ XX' + A' 2 + A + A' + A 2 + 3A 2 A /2 

-4AA ,2 -4A 2 A / ](A-1)(A'-1), (45) 



7 = -8(A-l) 2 (A'-l) 2 . (46) 

The leading order contribution to J^ a , in fact arises from the limiting region of x — > in 
the integral of Eq. (}4*3*|) . Therefore the higher order terms of 0(x 3 ) and 0(x 4 ) have been 
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ignored in the ^-function and the square root denominator in Eq. (|43|). The term f3x can 
be further neglected since it also is higher order smaller according to Eq. (|45j). Therefore, 
one has 



jA 



or 



= r (l, l) / dx 

9{x 2 - \{\ - l) 2 {\' - l) 2 ) 
4y / x 2 -±(A-l) 2 (A'-l) 2 



' -/-^d+A/A') 



(47) 



1 



^-|(A-1) 2 (A'-1)^ 



(48) 



Equation (j4*H|) can be evaluated as 



J^ = -V'(l,l)ln[(A-l)(A'-l)]. (49) 



Substituting the preceding result into Eq. (|22|) and carrying out the remaining integrations 
over z and z', one finally has 

Gt> = -^(1, 1)(P 2 " 4) 2 ln[(p 2 - 4)/24]- (50) 

In view of the fact that G^ a , in the preceding equation is totally similar to G^ a ? in Eqs. 
(JHHJ), and (|H4jl . it is curious that there seems no simpler way to derive it. 

Substituting the results for G% ly G^ a ,, and G® a , in Eqs. (JHIJ), (jHS)), and ((33) into Eq. 
(I26j) one finally obtains 

* -5(^)V-«r*--«« 

$> CTCT '(1, 1) + r CT '(l, -1) + !)]• (51) 

The contributions from the "direct" and "exchange" processes can be separately written as 



1 _ m% ln e 3 



r (D) 16vr 3 £; F 2£ F 

^T[2(W^(0)) 2 + (W aa ,(2k F )) 2 ], (52) 
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FIG. 2: Electron relaxation rate T in 2D 
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201 ] Dotted line: experimental data from 
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Ref. [5j]; dashed line: calculated one from Eq. (3) with RPA to the effective interaction W(q); 
solid line: calculated one from Eq. (3) with the effects of the local field and the finite well width 
on W(q) included. 



and 



T 



(ex) 



m 2 e P ln £ P 



16tt 3 E f 2E f 
[(^(0)) 2 + 2W aa (0)W aa (2k F ) 



(53) 



We empahsize that the above results are accurate only to the leading order of 0(£pln£ p ). 
In the case that the spin dependence of the effective interaction can be neglected, one has 
exactly the results shown in Eqs. (52) and (59) in I, respectively. The sum of and 
\/ T ( ex ) yields l/r e as given in Eq. (2) in the introduction. 



V. DISCUSSION AND SUMMARY 

The results in Eqs. (|15|) and (|16|) have been obtained with the approximation of the 
Thomas- Fermi screened Coulomb potential of Eq. (|14j) to the effective electron interaction. 
One can always resort to the more general expressions of Eqs. (fT2*j) and (fT3*j) if necessary. On 
the other side, in the 2D case, with which this paper is mainly concerned, no approximation 
has been made in the effective electron interaction W(q) except that it is assumed to be 
static. The local field effects and the finite well width effects on the effective interaction can 
be readily taken into account. We now estimate their contribution to the inverse quasiparticle 
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lifetime by using the form factor of F(q) = ~[1 + ^-(e qb — 1)] with b being the well 



qb<- qb 

width 



and the local field factor evaluated in Ref. [22]. These effects have been 
investigated earlier in Ref. . The clarification of the theoretical disagreement now enables 
us to definitely elucidate their contribution. Both of them are shown to yield in effect quite 
small corrections to the results calculated with random phase approximation (RPA) to W(q). 
This conclusion should hold in more general sense regardless of particular choice of the form 
factor and the local field factor, since both of them mainly effect the short-range behavior 
of the effective interaction while the inverse lifetime of a low energy quasiparticle is mainly 
determined by the long-range behavior of the effective interaction. The comparison with 
the experimental values of the electron relaxation rate T [ll[ [l^J llll Q from Ref. is 
illustrated in Fig. 2. The result calculated with the RPA to the effective interaction is also 



plotted in Fig. 2 for comparision. 



It seems that other factors must be taken into account in order to explain the difference 
between theory and experiment. One of the assumptions made in all previous work is that the 
couplings between electrons in different wells are weak and can be ignored. This assumption 
might require further justication for a barrier width being about 250A, while the width 
of each well being about 200A. Furthermore, higher order terms in electron interaction, 
usually not important at a density of r s ~ 1 where r s is the Wigner-Seitz radius, might not 
be simply ignored in this case, for they have been shown to contribute nontrivially higher- 
order lo gar ithmic factors. [3] Other factors which might also play a role have been metioned 
in Ref. [11]. Evidently, further theoretical effort is needed in order to fully understand the 
discrepancy between theory and experiment. 

In conclusion, we have calculated, in a rather different manner, the inelastic lifetime of an 
electron quasiparticle in an electron liquid. The results confirm those in Eq. (0), and Eqs. 
(0) and (JH1) obtained in our previous work, and consequently finally resolve the theoretical 
discrepancies in the literature. 
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